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ABSTRACT 

We solve numerically the ideal MHD equations with an external gravitational field in 
2D in order to study the effects of impulsively generated linear and non-linear Alfven 
waves into isolated solar arcades and coronal funnels. We analyze the region containing 
the interface between the photosphere and the corona. The main interest is to study 
the possibility that Alfven waves triggers the energy flux transfer toward the quiet 
solar corona and heat it, including the case that two consecutive waves can occur. 
We find that in the case of arcades, short or large, the transferred fluxes by Alfven 
waves are sufficient to heat the quiet corona only during a small lapse of time and in 
a certain region. In the case of funnels the threshold is achieved only when the wave 
is faster than 10 km/s, which is extremely high. We conclude from our analysis, that 
Alfven waves, even in the optimistic scenario of having two consecutive Alfven wave 
pulses, cannot transport enough energy as to heat the quiet corona. 
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1 INTRODUCTION 

Recent observations indicate the existence of Alfven waves 


Nakariakov et al. 1999) and torsional Alfven waves (Jess 


(Banerjee et al. 2007), kink waves (Aschwanden et al. 1999 


et al. 2009) in the solar atmosphere. These observations 


are based on data provided by current space missions, such 
as Transition Region and Coronal Explorer (TRACE), So¬ 
lar Dynamics Observatory (SDO), Solar Optical Telescope 
(SOT), X-ray Telescope (XRT) and Swedish Solar Telescope 


(SST) (jBanerjee et al. 2007||Tomczyk et al. 2007 

De Pontieu 

|et al. 2007 Okamoto et al. 2007| Jess et al. 2009 

). Moreover, 


it has been observed that the solar atmosphere is permeated 
by magnetic fields, which are organized in several structures, 
like small-scale magnetic flux tubes, loops, gravitationally 
stratified and magnetically confined structures (solar mag¬ 
netic arcades), weakly curved coronal magnetic flux tubes 
(coronal funnels), and others, and these structures may sup¬ 
port the propagation of differents kinds of MHD waves. 

We are interested in studying the role that Alfven waves 
play in the coronal heating problem, but first we mention 
the physical mechanisms that have been already proposed 
to study the problem. Theoretical models of coronal heat¬ 
ing mechanisms include Direct Current (DC) and Alternat- 
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ing Current (AC) models, which characterize the electrome¬ 
chanic coronal response to the photospheric driver that pro¬ 
vides the energy source for heating. The DC models involve 
Ohmic dissipation, magnetic reconnection, current cascad¬ 


ing and viscous turbulence (Sturrock & Uchida 1981 Parker 


1983 Milano et al. 1999 Gudiksen & Nordlund 2002), and 


the AC models involve wave heating by Alfvenic resonance 


(Hollweg 1991), resonant absorption (Davila 1987 Goossens 
et al. 199^ Ofman & Davila 1994 Erdelyi & Goossens 1994), 


phase mixing ( Parker 1991| De Moortel et al. 1999), cur¬ 
rent layers (Galsgaard & Nordlund 1996), MHD turbulence 
( Matthaeus et al. 1999[ Dmitruk et al. 2001) and cyclotron 


resonance (Hollweg 1986 Tu & Marsch 1997). In addition 


is it considered that some heating could also be produced 


by compressional waves (acoustic waves) or shocks (Kupe- 
rus et al. 1981). Another possible mechanism to heating the 


solar corona is by magneto-acoustic shock waves generated 
by chromospheric reconnection ( Sturrock 1999| Litvinenko 


1999 Sakai et al. 2000 Ryutova et al. 2001) 


The theoretical approach is centered in the analysis of 
numerical simulations modeling the propagation of Alfven 


waves in the solar atmosphere, e.g. in (Kudoh & Shibata 
1999) numerical simulations have been perfomed to study 


torsional Alfven waves propagation along an open magnetic 
flux tube in the solar atmosphere. The results suggest that 
the hot quiet corona, non-thermal broadening of lines, and 
spicules are all caused by the non-linear Alfven waves. In 
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(Del Zanna et al. 2005), it has been studied numerically the 


propagation and evolution of Alfvenic pulses in the solar 
coronal arcades. The results show a propagation of the ini¬ 
tially localized pulses within the solar arcade, due to the 
variations in the Alfven speed with height, and therefore 
an efficient damping of the amplitude of the oscillations. In 


(Murawski & Musielak 2010) small-amplitude Alfven waves 


are modeled by solving the ID wave equation and it was 
found that as a result of cut-off frequencies, these waves are 
not able to propagate freely into the upper regions of the so¬ 
lar atmosphere, undergoing reflection towards lower layers. 


On the other hand, in (Ebadi et al. 2012 ) the dissipation 


of Alfven waves due to phase mixing in a stratified environ¬ 
ment of solar spicules has been presented . The results show 
that in the presence of stratification due to gravity, damping 


takes place. In (Chmielewski et al. 2013), it has been inves¬ 


tigated the role of impulsively generated non-linear Alfven 
waves in the observed non-thermal broadening of some spec¬ 
tral lines in solar coronal holes. The results were that the 
non-linear Alfven waves with amplitude Av = 50 km/s are 
the most likely candidates for the non-thermal broadening 


of line profiles in the polar coronal hole. In (Chmielewski 


et al. 2014a) a 2.5 D numerical model is used to study the 


propagation and reflection of impulsively generated Alfven 
waves within a solar magnetic arcade; the results show that 
the propagation of Alfven waves is affected by the spatial de¬ 
pendence of the Alfven speed, which leads to phase mixing 
that is stronger for more curved and larger magnetic arcades. 


In (Chmielewski et al. 2014b) we found a study of the propa¬ 


gation of impulsively generated linear and non-linear Alfven 
waves in weakly curved coronal magnetic flux-tubes, or coro¬ 
nal funnels, and it was found that non-linear Alfven waves 
may carry enough energy as to heat the coronal funnels and 
power the solar wind that originates in these funnels. 


Within all this rich and well investigated scenario we 
raise again the question of whether Alfven waves are capable 
of transferring enough energy as to heat the quiet corona. 
Most of the aforementioned studies focused on this question 
are based on numerical simulations of a single pulse. In this 
paper we simulate the propagation of two pulses in order to 
study the possible effects of the interference in the process 
of energy transfer to the corona. 


We study the possible effects of interaction between two 
impulsively generated Alfven waves in two configurations of 
magnetic field: i) solar magnetic arcades, that we assume to 
be isolated and fixed in time so that they may apply to hold 
on a quiet corona, and ii) coronal funnels. We are interested 
in the implications of Alfven waves in the heat transfer into 
the corona. We simulate the interaction by launching two 
pulses separated by a given interval of time between 0 and 50 
s, considering that otherwise the two pulses are considered 
to be independent. 

The paper is organized as follows, in section we de¬ 
scribe the model of solar atmosphere, the numerical meth¬ 
ods we use, and the way we measure the energy flux carried 
by Alfven waves. In section we present the results of the 
numerical simulations for various experiments. Finally, in 
section we present the conclusions. 


2 MODEL AND NUMERICAL METHODS 

In this section we describe the equations modeling the 
plasma, the photosphere-corona interface model, the dimen¬ 
sionless equations solved by our code, the initial conditions 
of the Alfven wave pulses and the numerical methods used 
to solve the evolution equations of the system. 


2.1 The system of equations 

We consider a gravitationally stratified solar atmosphere, 
which is described by a plasma obeying the ideal MHD sub¬ 
mitted to an external constant gravitational field. We choose 
to write down these equations in a conservative form, which 
is optimal for our numerical methods: 


| + V^(,v)=0, 

dipw) „ / B (g) B ^ . 

^ ^ + V • (pv (g) V-h Ipt) = pg, 


dt 


po 


ri'F' R 

^ + V • ((S + pt)v-(v B)) = pv g, 

Ot Po 

^B 

-^ + V-(v(g)B-B®v) = 0, 

V-B = 0, 


( 1 ) 

( 2 ) 

(3) 

(4) 

(5) 


1^2 

where p is the plasma mass density, pt = p + is the to¬ 
tal (thermal + magnetic) pressure, v represents the plasma 
velocity field, B is the magnetic field, po is the vacuum per¬ 
meability, I is the unit matrix and E is the total energy 
density, that is, the sum of the internal, kinetic, and mag¬ 
netic energy density 


E = 


P . 

7-1 2 2po ■ 


( 6 ) 


For the fluid we consider an adiabatic index 7 = 5/3. The 
system of equations 0-0 is closed with the ideal gas law 


kB rr. 

p = —pT, 

m 


(7) 


where T is the temperature of the plasma, m = firriH is 
the particle mass, which is specified by a mean molecular 


weight value of p =1.24 for a fully ionized gas (Chmielewski 


et al. 2014a) and uih is hydrogen’s mass, ks is Boltzmann’s 


constant. The gravitational source terms on the right hand 
side of equations 0 and 0 is given by g = [0, —p, 0] with 
magnitude g = 274 m/s^, which holds near the photosphere. 


2.2 Dimensionless MHD equations 

In order to carry out our numerical calculations we rescale 
the system ( T|[5 ) according to the following conventions. We 
choose the length scale /o, plasma density scale po, and mag¬ 
netic field scale Bq to fix the unit of time in terms of the 
Alfven speed 


_ _ ^0 — ^0 

Vo = VA,0 = , ^ to = -. 

^//M)P0 Vo 


( 8 ) 
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Table 1. Units used in this paper 


Quantity 

cgs units 

lo 

10® cm 

VO 

10® cm-s“^ 

to 

1 s 

po 

10-15 gr.cm“® 

Bo 

14.472 G 


a non-trivial magnetic field ^-component Bz that vary with 
X and y are allowed. This kind of model is called a 2.5 dimen¬ 
sional (2.5D) model ( [Murawski Musielak 2010[ Murawski 


^^aqarashvil^^2010j^Miirawsk^^e^^l^201lUWoloszkiew^ 

et al. 2014| |Chmielewski et al. 2014a| |Chmielewski et al. 
2Q14b| ). 


Given the solar atmosphere is in static equilibrium, it 
should have a force-free and current-free magnetic field given 
by the condition 


Using these basic scale factors /o, ^o, to we define dimen¬ 
sionless independent variables and associated differential op¬ 
erators as follows 



V = /oV, t = - 
to 


— = / — 
dt 


(9) 


and consequently the corresponding dimensionless state 
variables are defined by 


V X Be = 0 ^ (V X Be) X Be = 0, (18) 

where the subscript e stands for equilibrium quantities. The 
choice of the magnetic field satisfying (18) is general and 
depends on the physical problem under study. As a result 
of applying the static equilibrium and the magnetic current- 
free conditions to equation it- the pressure gradient is bal¬ 
anced by gravity: 


P = 


Po 


Vo 


P = 


povl' 


E = 


povl' 




lo 


( 10 ) 


Finally the dimensionless equations ruling our system are 


| + V.(/,v) = 0, 


d{p\ 


dt 


+ V • (pv (g) V - B (g) B + Ipt) = pg, 


f)E - - - - 

^ + V-((S + pt)v-B(v-B))=pv-g, 

^B - - - 

-^ + V-(v®B-B®v)=0, 

V - B = 0. 


( 11 ) 

( 12 ) 

(13) 


The dimensionless version of the equation of state 0 
obtained using the definitions of p and p 


(14) 

(15) 
is 


-Vpe + Peg = 0. (19) 

In addition, following ([Murawski Musielak 201(5) Mu-| 


rawsk^^^^agarashv^ ^ 010| [Murawski et al. 2011 
Woloszkiewicz et al. 26l4| fChmielewski et al. 2014a 


Chmielewski et al. 2014b) we use the ideal gas law (0 and 


the fact that hydrostatic equilibrium will only hold along 
the y direction, so that the equilibrium gas pressure takes 
the form: 


Pe{y) ^Prefe^P 
and the mass density 



( 20 ) 


where 


Pe{y) 


Pejy) 

gHy) ’ 


( 21 ) 


ks rr 

P = —pT 
m 


p 


Pb 

v^m 


pT. 


(16) 


and the dimensionless temperature T and equation of state 
are 


T = T 


Pb 


P = pT. 


(17) 


In this way, the ideal MHD system is attached to 

the scale length, magnetic field, density and time factors /o, 
Bo, Po and to - In our simulations we choose the scale factors 
in cgs units as specified in Table 


AO , teKM, (22) 

mg 

is the pressure scale-height, and pref represents the gas pres¬ 
sure at a given reference level, that we choose to be located 
at yref — 10 Mm, because in the temperature model above, 
this is the location at which the corona region starts. 


2.3.1 Model for isolated solar magnetie areades 


We model an isolated solar magnetic arcade with a magnetic 
flux-tube model originally developed by (Low 1985) in three 
dimensions. This model considers a magnetic flux function 


2.3 Model of the static solar atmosphere 

We specifically center the domain of analysis as one covering 
part of both the chromosphere and the corona. We consider 
the atmosphere to be static {d/dt = 0) and study the evolu¬ 
tion on a finite xy domain, where the x is a horizontal coordi¬ 
nate and y labels height. All the state variables depend on x 
and y, but a non-zero component of the velocity field Vz and 


where 6 is a constant that determines the vertical location 
of a magnetic pole, which is set to b = —5 Mm and Bref 
is the magnetic field strength at yref - The magnetic field is 
thus 
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Figure 1. (Top) Magnetic field lines for an open arcade with 
parameters Bref = 14.472 G, b = —5 Mm and i/ref = 10 Mm. 
(Bottom) The associated Alfven speed in Mm/s. 


Be(a;,j/) = V X (A(x,y)z) = 


dA{x,y) ^ _ dA{x,y) , 
dy ^ dx 


(24) 


where specifically the equilibrium magnetic field components 
Bex(x,y) and Bey(x,y) are 


Ae(x, y) = Kb Bref sin ( 


exp 


y yref 
Kb 


(28) 


where Bref is again the magnetic field at the reference level 
yref = 10 Mm. Unlike (Chmielewski et al. 2014b) we use 
a sin instead of a cosine in order to be consistent with our 
convention of coordinates. Here Kb = 2L/7r denotes the 
magnetic scale-height and L is a half of the magnetic arcade 
width. In this case we want to model weakly expanding coro¬ 


nal funnels, then we fixed L = 75 Mm (Chmielewski et al. 


2014b). For this choice, the magnetic field lines are weakly 


curved and represent the open and expanding field lines sim¬ 
ilar to coronal holes. Finally, the magnetic field reads 


Be(a;,j/) = V X {A{x,y)z) = 


a^(a:,2/). dA{x,y)^^_ 

dy dx 


and the equilibrium components Bex{x, y) and Bey{x, y) are 
explicitly 


Bex {x, y) 


Bey{x,y) 


R • I ® 

— Bref Sm I - 

^Bref eOS 


\Kb 

fe) 


exp 

exp 


y yref 
Kb 

_ y ~ yref 
Kb 


(30) 

(31) 


On the other hand, in this model the Alfven speed varies 
only with y because the density does: 


ca( 2 /) = 


■Bfa; + Bly 
nopeiy) 


y-Vref 

Bref^ 

VPoPeiy) 


and the equilibrium sound speed is 


(32) 


Bex(x,y) = - 

Beyix,y) = - 

The value of Bref 
speeds 


2x{y - b){b - yreff „ 

(a:2 - (b - j/)2)2 
{x'^-\-{b-yf){b-yreff 


B, 


ref- 


(25) 


(26) 


(y^ — 2by — + 5^)^ 

is choosen such that the Alfven and sound 


CA{x,y) = 


Blx + Bly 
poPeiy) 


Ay) = 


_ hpeiy) 


Pe{y) ’ 


satisfy the constraint CA{0,yref) = 10cs(yre/), because in 
the solar corona the magnetic pressure is bigger than the 
fluid pressure. This constraint gives a value of Bref = 14.472 
G. The magnetic field components resulting from equations 
(25) and (26) are displayed in Fig. (top). Such magnetic 
configuration corresponds to an isolated asymmetric mag¬ 


netic arcade as proposed in (Chmielewski et al. 2014a). As 


described in equation (27) Alfven speed depends on x and 


y, which means that it is non-isotropic as is also shown in 
Fig. a (bottom). This effect will influence the dynamics of 
the wave. 


2.3.2 Model for coronal funnels 


Following (Chmielewski et al. 2014b), we assume the mag¬ 
netic flux function has only a z component and the following 
expression 


cs{y) 



(33) 


In addition, Bref must satisfy the condition CAiyref) — 
10cs(yre/), and once again we set Bref = 14.472 G as in the 
case of magnetic arcades. The resulting magnetic field and 
Alfven speed are shown in Fig. Due to the symmetries, 
unlike the arcade, in this case the Alfven speed depends only 
on y. 


2.4 Plasma temperature model 


The temperature field is assumed to obey the C7 model 


(Avrett & Loeser 2008). This is a semiempirical model of the 


chromosphere, with the temperature distribution adjusted 
to obtain optimum agreement between calculated and ob¬ 
served continuum intensities, line intensities, and line pro¬ 


files of the SUMER (Curdt et al. 1999) atlas of the extreme 


ultraviolet spectrum. 

We set the temperature profile by interpolation of the 
data in ( Avrett fc Loeser 2008| into our numerical domain, 
and the result is shown in Fig. The main properties of 
this profile are that T varies nearly three orders of magni¬ 
tude across the transition between the chromosphere an the 
corona, ranging from 5000 K at y = 1.5 Mm up to 10® K at 
y — Mm. For bigger y the temperature profile is assumed 
to be constant. 
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x(Mm) 



Figure 2. (Top) Magnetic field lines for coronal funnels with 
parameters = 14.432 G, Href = 10 cind = ISO/tt Mm. 

(Bottom) Alfven speed as a function of y. 


The pressure pe (y) is integrated numerically using equa¬ 
tion (20) and Pe{y) is calculated using equation (21). The 
two profiles are shown in Fig.^ where the important gradi¬ 
ents at the transition region can be seen. 


2.5 Initial conditions 


Using the previous experience in ( [Murawski Musielak 
2010| [Chmielewski et al. 2Q13| |Woloszkiewicz et al. 2Q14| 


Chmielewski et al. 2014a[ [Chmielewski et al. 2014b ) the 
Alfven waves are considered to be perturbations of the 
2 :—component of the velocity field Vz with a given profile. 
In our cases, both for the arcades and the funnels we use 


the velocity profile in (Chmielewski et al. 2014a) consider¬ 


ing a Gaussian profile of the form 



Figure 3. Equilibrium temperature profile for the C7 model 
interpolated into our numerical domain. 




_ (x-xq)^ _ (y-yp)^ 

v,{x,y,t = 0) = A,e , (34) 

where Ay is the amplitude of the pulse, (xo, yo) is its initial 
position, Wx and Wy denote the width of the pulse along x 
and y. 

In this paper we consider the interaction between two 
pulses, one launched at initial time and a second pulse 
launched after a given delay lapse r. We study the effects 


Figure 4. Equilibrium gas pressure (top) and mass density (bot¬ 
tom) for the C7 model. 


of the value of the delay lapse on the dispersion of the to¬ 
tal energy of the pulses and their capability of transporting 
energy toward the corona. 
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2.6 Numerical methods 


The equations are solved numerically using our 


own numerical code, which is the Newtonian version of 


CAFE (Lora et ah 2014), that solves the MHD equations on 


a uniform cell centered mesh in unigrid mode, and has been 
applied in relativistic scenarios including the accretion on 
black holes (Lora et ah 2014b Lora-Clavijo & Guzman 2013 


Cruz-Osorio, et ah 2012 Guzman & Lora-Clavijo 2011). The 


integration in time uses the method of lines with a third or¬ 


der Runge-Kutta time integrator (Shu & Osher 1989). The 


right hand sides of the MHD equations are discretized using 
a finite volume approximation, together with a High Resolu¬ 
tion Shock Capturing method. In all the simulations in this 


paper, this method uses the HLLE (Harten P. et al. 1983) 


approximate Riemann solver formula in combination with 
the second order linear piecewise reconstructor MINMOD. 

The numerical evolution of initial data involving 
Maxwell equations leads to the violation of the divergence 
free constraint equation § , developing as a consequence un¬ 
physical results like the presence of a magnetic net charge. 
Among the several methods available to control the growth 


of the constraint violation (Toth G. 2000), in our code we 


use a version of the constrained transport method (CT) pro¬ 
posed in (Evans & Hawley 1988), which is based on the use 
of the fluxes computed with the conservation scheme itself. 


This algorithm is known as flux-CT (Balsara 2001), which 
maintains the constraint small. 

We use a 2D domain, and the boundary conditions con¬ 
sist in fixing all the evolution variables at all four boundaries 
to their equilibrium initial values. 


2.7 Alfven waves and coronal heating 

In order to estimate the influence of Alfven waves as carriers 
of energy to the solar corona, which eventually would pro¬ 
duce the heating of the quiet corona or the acceleration of 
the solar wind, it is necessary to have flux estimates toward 
the corona. One way to measure the flux energy carried by 
Alfven waves is to consider the Poynting flux given by 


Fp = —B X (v X B) 
Mo 


(35) 


from which we only take the y—component that carries 
energy toward the upper shells. Another usual measurement 
involves the WKB approximation, which approximates the 


& Shibata 1999): 


total flux formula as (Browning 1991 Hollweg 1991 Kudoh 


' P '^zCA, 


(36) 


where p is the density, ca is the local speed of the Alfven 
waves and Vz is the Alfven speed, measured at a particular 
height. Here we measure both fluxed in order to compare. 

A way to estimate whether the Alfven waves contribute 
or not to heat the quiet corona is using the threshold in 
(Withbroe & Noyes 1977), indicating that when the flux 


exceeds the value ~ 3x10^ erg/s/cm^, the contribution to 
heating is considerable. Usually what is done is that Fwkb 


is measured at a specific point (Kudoh & Shibata 1999 


Chmielewski et al. 2014b). Here we measure both fluxes not 


at a single point, but in order to have a bigger picture we 
scan these values in three different points specified for each 
case. 


3 RESULTS OF NUMERICAL SIMULATIONS 

We present the results of the energy flux transfered to the 
corona due to the interaction of two waves evolving on ar¬ 
cades and on funnels. In all cases we estimate the effects 
of the time elapsed between the two pulses r. For this we 
only consider impulsively generated Alfven waves because 
their higher amplitude corresponds to a better transport of 
energy. 


3.1 Interaction between two Alfven waves on a 
solar magnetic arcade 

We set the simulation domain to [1,29] x [1,15] Mm^, covered 
with 1200x600 cells. We consider large and short arcades, 
and the length of the arcades depends on the curvature of the 
magnetic field lines: a) For the large arcades we use xo = 2 
Mm, and the Alfven waves propagate higher up along longer 
and more curved magnetic field lines, b) for short arcades 
we set xo = 3 Mm, and the inclination of the magnetic field 
lines is much larger and the magnetic field lines are less 
curved. 


3.1.1 Interaction of two pulses in a large arcade 


In this case the amplitude of the initial Gaussian velocity 
pulse is = 3 km/s, the widths are Wx — 0.2 Mm, Wy — 0.1 
Mm and is located initially at (xo,Mo) =(2 Mm, 1.75 Mm) 
in ( [34| . We show the results for the particular case in which 
the second pulse is launched 30 s after the first pulse in Fig. 
??. Initially the first pulse decouples into two pulses mov¬ 
ing in opposite directions. One pulse propagates downwards 
very fast, and the second propagates upwards through the 
chromosphere and transition region to the solar corona, the 
part of the pulse that propagates upwards reaches the tran¬ 
sition region and accelerates due to an increase of the local 
Alfven speed as expected for the case of a single pulse in 


(Chmielewski et al. 2014a). We launch the second pulse at 


time t = 30 s, this pulse clearly affects the propagation of 
the first pulse. The interaction between the two pulses pro¬ 
duces reflections as shown at time t = 76 s. Finally by time 
t = 106 s, the waves reach a maximun distance x where they 
are reflected back to the solar surface. This process shows 
that phase mixing and reflections grow due to the interac¬ 
tion. Finally, in order to verify our evolution we show at 
the bottom of Fig. ?? the violation of the magnetic field 
divergence free constraint. 

We perform a series of simulations like this one, for 
various time lapses between the first and the second pulses 
T. We measure Fpy and Fwkb at three different points and 
show the measurements in Fig. We plot a line indicating 
the threshold 3 x 10^ erg/s/cm . The results indicate that 
the Poynting flux measured at the detector located on the 
left overtakes the threshold during a small lapse of time. The 
reason is that at this detector the field is more intense. This 
cannot be observed at other detectors and according to the 
Fwkb the threshold is never achieved. 
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Figure 6. Energy fluxes Fpy (top) and Fwkb (bottom) for the propagation of two pulses in a large arcade using three detectors 
(from let to right) located at (7.25, 7) Mm, (14.5, 7) Mm and (14.5,10) Mm, for r = 0, 2, 5,10, 20, 50 s. The continuous line indicates the 
threshold over which the flux contributes to the coronal heating. 


3.1.2 Interaction in a short arcade 

In a similar way, we study the propagation of the same ve¬ 
locity Gaussian pulses with Ay = 3 km/s, Wx = 0.2 Mm, 
Wy = 0.1 Mm and (xo,yo) =(3 Mm, 1.75 Mm) on a short 
arcade. In Fig. ?? we show a particular case with r = 20 s. 
Again, the initial pulse decouples into two pulses moving in 
opposite directions. At time t = 16 s we only show the part 
of the pulse that moves upwards without interaction. This 
pulse is affected by the transition region and phase mixing, 
which produces the appearance of a fanning with negative 
amplitude. We launch the second pulse at time t = 20 s, 
which produces interference at time t = 32 s. By t = 40 s 
the second pulse suffers the effect of the density decay in 
the transition region, and the fanning becomes elongated. 
At time t = 52 s the waves reach a maximum distance in 
w, where they suffer strong reflections. Part of these waves 
return to the solar surface, so as the energy. 

We perform the same experiments as in the previous 
case for various time lapses between the first and the sec¬ 
ond pulses. Again we use three detectors to measure Fpy 
and Fwkb shown in Fig.[^ Like in the previous case, only 
the Poynting flux reports values above the threshold in the 
detector at the left. 


3.2 Interaction between two Alfven waves in a 
coronal funnel 

In this case we set the simulation domain to be [-5,5] x[l, 
81] Mm^, covered with 300x2400 cells. Like in the cases 
of arcades, we launch two velocity pulses with profile 


separated by an interval of time r on the funnel magnetic 
field configuration defined in |2.3.2[ 


3.2.1 Interaction between two Alfven waves with 
amplitude Ay — It) km/s 

The typical behavior of the pulses is shown in Fig. for 
T = 10 s and pulses with parameters = 10 km/s, Wx — 
Wy — 0.2 Mm, {xq =0, yo =1.75) Mm. Initially the first 


pulse decomposes into two propagating waves (Chmielewski 


et al. 2014b), one pulse that moves upwards and a second 


pulse that moves downwards. The amplitude of the pulse 
that moves downwards decays with time and practically dis¬ 
appears. We show the first pulse at time t = 5 s, this pulse 
experiences an acceleration at the transition region due to 
an increment of the local Alfven speed ca • The shape of the 
pulse becomes ellipsodial elongated along the vertical direc¬ 
tion. Then the second Gaussian pulse appears at time t = 10 
s. Still by t = 15 s the effect of the second pulse is not signif¬ 
icant to the shape of the first pulse. However, at that time 
the pulse that moves upwards suffers another elongation. 
At later times in the evolution, we can see reflections and 
phase mixing because of the effect of the transition region 
in the pulses that move upwards at time t = 30 s. The two 
pulses have already reached the solar corona. At time t = 60 
s, a signal of Vz with negative amplitude reaches the solar 
corona, and the pulses that move upwards continue moving 
vertically. 

The fluxes measured at various detectors are shown in 
Fig. |10| In this case the Poynting flux never reaches the 
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Figure 8. Energy fluxes Fpy and Fwkb for the propagation of two pulses in a large arcade using three detectors at (7.25,4.5) Mm, 
(14.5, 4.5) Mm and (14.5, 7) Mm, for r = 0, 2, 5,10, 20, 50 s. Again we show the threshold. 


threshold for any of the combination of pulses whereas, un¬ 
like the arcade cases, the Fwkb flux exceeds the threshold 
during a small lapse of time for the case r = 5s, measured 
by a detector located along the central vertical line (0,10). 


3.2.2 Interaction between two Alfven waves with 
amplitude Ay = 40 km/s 

In order to show the contribution using one of the highest 
amplitudes found in literature, we prepare a case of an im¬ 
pulsively generated and high amplitude pulse with parame¬ 
ters Ay = 40 km/s, Wx — Wy = 0.2 Mm, (xo =0, yo =1.75) 
Mm (Chmielewski et al. 2013 Chmielewski et al. 2014b). 
We show snapshots of the evolution of the wave in Fig. m 

In Fig. we show Fpy and Fwkb ut various detec¬ 
tors. We And that the two detectors located at the center of 
the domain measure fluxes, both Fpy and Fwkb over the 
threshold for various values of r in a generic way. 

As mentioned before, this is one of the highest value 
used to model the propagation of Alfven waves, however it is 
not clear from observations that such high speeds have been 
observed and associated to funnel type of configurations so 
far. Therefore we have only shown that in the case that such 
high speeds are observed, they may actually contribute to 
heating the corona. 


corona, but two possible pulses that may be launched very 
close in time. 

We found that impulsively generated linear Alfven 
waves in the solar magnetic arcades, with slow pulses, in 
large or short arcades, carry enough energy as to heat the 
solar quiet corona only during a small lapse of time for a 
very particular combination of pulses, which indicates the 
threshold is not overtaken for generic combination of pulses. 

In the case of the propagation of impulsively generated 
non-linear Alfven waves in the coronal funnels, the waves 
with amplitud = 10 km/s do not carry enough energy 
flux as to heat the solar quiet corona, even if two consec¬ 
utive pulses are launched. On the other hand, the energy 
transferred overtakes the threshold by orders of magnitude 
in the case of waves with an amplitude of Ay = 40 km/s, 
which we cannot be sure that exist but has been mentioned 
in existing literature. 

We conclude from our analysis, that impulsively gen¬ 
erated Alfven waves, even in this optimistic scenario, are 
unlikely candidates to heat the quiet corona. 


4 CONCLUSIONS 

We have proposed an optimistic scenario where not only a 
single Alfven wave pulse can contribute to heating the quiet 
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Figure 9. Spatial profiles of the interaction of two non-linear pulses 
times t = 5 s, t = 15 s, t = 30 s, t = 60 s. In this case r = 10 s. 
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Figure 11. Spatial profile of the interaction of two non-linear pulses Vz{x,y,t) for Ay = 40 km/s, Wx = Wy = 0.2 Mm displayed at 
times t = 5 s, t = 15 s, t = 30 s, t = 60 s. In this case r = 10 s. 






































